Single-particle excitations in a trapped gas of Fermi atoms in the 

BCS-BEC crossover region 

Y. Ohashi 1 * 2 and A. Griffin 2 

1 Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan, 
2 Department of Physics, University of Toronto, 
Toronto, Ontario, Canada M5S 1A7 
(Dated: February 2, 2008) 



1 



Abstract 

We investigate the single-particle properties at T = of a trapped superfluid gas of Fermi atoms 
with a Feshbach resonance. A tunable pairing interaction associated with the Feshbach resonance 
leads to the BCS-BEC crossover, where the character of superfluidity continuously changes from 
the BCS-type to a BEC of composite bosons (consisting of a superposition of a condensate of 
Cooper-pairs and molecular bosons). In this paper, we extend our previous work for a uniform 
superfluid Fermi gas [Y. Ohashi and A. Griffin, Phys. Rev. A 67, 063612 (2003)] to include the 
effect of a harmonic trap. We do not use the local density approximation (LDA), but directly 
solve the Bogoliubov-de Gennes coupled equations. Using these equations, we find self-consistent 
values for the spatially-dependent local density n(r) as well as the composite BCS order parameter 
A(r), the latter describing both the Cooper-pair and molecular condensate contributions. Using 
these results, we calculate the single-particle density of states in the crossover region, and from this 
determine the true single-particle energy gap (E g ) of the trapped Fermi superfluid at T = 0. This 
is associated with the in-gap (or Andreev) states in the low density region at the edge of the trap. 
We calculate the laser-induced current into another hyperfine state, as measured in recent 
rf-spectroscopy experiments. This rf-spectrum gives a direct probe of the quasiparticle spectrum. 
We show how the high-energy part of I(u>) gives information about A(r = 0) at the center of the 
trap (which is comparable to the Fermi energy ep m the crossover region) . More generally, we show 
that I(uj) is very dependent on the spatial profile of the pair potential A(r) which is used. We also 
emphasize that the narrow "unpaired atom" peak in the rf-data gives information about E g and 
the low-energy (^C ep) in-gap states of a Fermi superfluid. While our calculations are limited at 
T = 0, we use them to discuss the recent Innsbruck data and the LDA calculations of Torma and 
co-workers. The LDA, while useful, can lead to an incorrect physical picture of the low density 
surface region of the Fermi superfluid. 

PACS numbers: 03.75.Ss, 03.75.Kk, 03.70.+k 
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I. INTRODUCTION 



Recently, the BCS-BEC crossover in a trapped gas of Fermi atoms near a Feshbach 



n a trappea 

mw 



Molecular bosons associated with a 



resonance has attracted much attention 
Feshbach resonance can mediate a tunable pairing interaction between atoms, which becomes 
stronger with decreasing threshold energy (denoted by 2v) of the Feshbach resonance^, 
The BCS-BEC crossover can be studied as a function of 2z/, with the character of 
superfluidity continuously changing from the conventional BCS-type of Cooper-pairs to a 
BEC of composite bosons (consisting of a superposition of Cooper-pairs and molecules), as 
one approaches the strong-coupling regime jl|, 0, 0, Q|- This situation contrasts with the 
"classic" crossover physics originally studied in the superconductivity literature 

HQ 



13, 



14, 



15(, where Cooper-pair bosons are always dominant in the whole crossover 



regime, from the BCS to the BEC limits. Recent experiments on ultracold gases of 40 K 
and 6 Li have used this tunable interaction near the resonance ^|, to produce a large 



number of 
length) 



Dound states or molecules when a 2b > (a 2b is the two-body s-wave scattering 



19j. Using this tunable interaction, a BEC of these molecular bosons has been 



observed in K and Li 



2l|, and more recently, evidence for superfluidity was found in 



the BCS side of the crossover regime (loosely defined as a? b > 0) 

nn 
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rnuictity was ton 



Both 



single-particle excitations [27|] and collective modes[2q, |28( have been experimentally studied 
in the crossover regime in the case of 6 Li. 

In this paper, we present a detailed study of BCS-type single-particle excitations in the 
BCS-BEC crossover of a trapped gas of Fermi atoms with a Feshbach resonance. This 
extends our previous work for a uniform gas 3] in a major way, since we now include the 
effect of discrete eigenstates due to confinement in a harmonic trap. However, in contrast 
to Ref. 0| , we limit our discussion to T = in this paper. Because of the inhomogeneity of 
a gas due to the trap potential, the single-particle threshold excitation gap E g is no longer 
simply related to the position-dependent superfluid pair potential A(r). This contrasts with 
the case of a uniform Fermi superfluid 0, Is, Q], where E g and the order parameter A are 
related in a very direct way (E g = |A| for jj > 0, E g = \J fi 2 + |A| 2 for /i < 0, where /i 
is the chemical potential of the fermions). We compute A(r) self-consistently by solving 
the Bogoliubov-de Gennes (BdG) coupled equations and then use it to calculate the single- 
particle density of states N(u>) for a trapped superfluid gas. We present results through the 
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BCS-BEC crossover by varying the threshold energy 2v. We also calculate the rf-tunneling 
current, which can be used to probe the spectrum of the single-particle excitations. We 
compare our T = results with recent experimental data on 6 Li 27[. This data agrees with 



the calculations of Torma and co-workers 



29|, using a local density approximation (LDA). 



While the LDA is useful, we feel that it gives an incorrect picture of the surface region of 
the trapped superffuid. 

This paper is organized as follows: In Sec. II, we introduce the usual coupled fermion- 
boson Hamiltonian including an isotropic harmonic potential and derive the mean-filed BdG 
equations. The single-particle Green's functions in a trap are expressed in terms of the 
solutions of the BdG equations in Sec. III. In Sec. IV, we discuss the BCS-BEC crossover 
in the cases of broad and narrow Feshbach resonances. We review different definitions of 
the s-wave scattering length. In Sec. V, we calculate the equilibrium properties in a self- 
consistent way, such as the chemical potential /i (which plays a crucial role), the atomic 
density profile n(r) and spatial variation of the composite order parameter A(r). The 
single-particle density of states N(u) for different values of 2v is calculated and discussed 
in Sec. VI. In Sec. VII, we compare our results (both BdG and LDA) with the rf-tunneling 
spectroscopy data obtained in recent experiments on 6 Li. Some of our results were briefly 
reported earlier in Refs-I^. I3I | . 



II. EXTENSION OF THE BOGOLIUBOV-DE GENNES EQUATIONS (T = 0) 

We consider a two-component Fermi gas with a Feshbach resonance trapped in a harmonic 
potential, using the coupled fermion-boson model 

H = £ / dr*l(T)[~- f i + V* 9 (T)]* <r (T)-U [ »{(r)*{(r)^(r)* t (r) 



dr$ f ( 



2M 



+ 2u - 2// + <ap(r)l$(r) + g r / dr U^r)*^)*^) + h.c 



(2.1) 



Here ^(r) is the fermion field operator with pseudo-spin a =f , |. The Bose quantum field 
operator $(r) describes molecular bosons associated with the Feshbach resonance. U is a 
non-resonant interaction, which we take attractive {—U < 0). The Feshbach resonance g r 
describes resonance between one molecule and two Fermi atoms. The effect of this resonance 
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is controlled by adjusting the energy 2v of the molecules, also referred to as the threshold 
energy. (To avoid confusion, we note that in many recent papers, this threshold energy is 
denoted by v.) Since a molecule consists of two Fermi atoms, we take M = 2m and impose 
the conservation of the total number of atoms, 

N = Np + 2N M 

= J drnp(r) + 2 / drnu(r) 

= W <*r(*t (r)* CT (r)) + 2 / dr($t(r)$(r)>. (2.2) 

This constraint has already been taken into account in Eq. (j2.1j) . with the Fermi chemical 
potential fi and Bose chemical potential „„ s 2/iQ. <, p (r) and V« (r) are the harmonic 
trap potentials for Fermi atoms and Bose molecules, respectively, which are assumed to be 
isotropic 



VL P = ^lr\ Cap=o^oV 2 - (2-3) 



In addition, in this paper, we assume that the atoms and molecules feel the same trap 
frequency u = ujqm (which correctly describes recent experiments). 

To study. the BCS-BEC crossover phenomenon, we extend the theory at T = developed 
by Leggettp] in the context of superconductivity to include a Feshbach resonance as well as 
the effect of a harmonic trap potential. The key point of this theory is to solve the mean-field 
gap equation for the order parameter together with the equation for the number of particles, 
which determines the Fermi chemical potential \x. In the BCS-BEC crossover, /i decreases 
from the usual BCS limit given by the Fermi energy e-p and can become negative. The 
thermal fluctuations in the Cooper-channel and the thermal excitations of Bose condensate 
fluctuations, crucial in considering the BEC-BEC crossover at finite temperatures close to 



T, 



10 



12 



13 



14 



151 ] - are not important at T = 0. Such fluctuations will not be 
considered in the present paper. 

The gap equation is obtained from the mean-field approximation for Eq. ()2.1j) in terms 
of the BCS Cooper-pair condensate A(r) = f/(\I / j(r)\I/^(r)) as well as the molecular BEC 
condensate 4>m (r) = (<&(r)). The Hartree-Fock-Bogoliubov (HFB) mean- field Hamiltonian 
for the fermions is given by 

V 2 U 



#hfb = E / «t« [-^ - 2-tifOO + V£p(r) - ft 
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^{(r)^{(r) + h.c. 



+ g r I dr<j) M (r) [tfj(r)*I(r) + h.c. 



rt, 



E / »t(r) [-^ - |n F (r) + V£ ap (r) - /i 



drA(r) ^{(r)^{(r) + h.c. 



(2.4) 



In Eq. (|2.4j) . the Cooper-pair order parameter A(r) and the molecular condensate 0m ( r ) 
are taken to be real, without loss of generality. np(r) = ^2 a (W a (r)^ a (r))) is the local 
number density of Fermi atoms. Equation (|2.4jl clearly has the same form as the usual 
BCS-Gor'kov Hamiltonian, but now with the Cooper-pair order parameter A(r) replaced 
with the composite order parameter, 



A(r) = A(r) -<? r </> M (r). 



(2.5) 



Because of the spherical symmetry of our harmonic trap, A(r), <pM{ r ), A(r), and rip(r) only 
depend on |r|. 

Since we are only discussing T = in this paper, the molecules described by $(r) are 
all Bose-condensed. As noted above, we ignore excitations of this molecular condensate. 
Thus the HFB mean-field Hamiltonian in Eq. (J2.4)) does not involve the dynamics of the 
molecular condensate. The latter only enters through the equilibrium value A(r). 

Superfluidity in our coupled fermion-boson model is characterized by the two broken- 
symmetry order parameters, i.e., the BCS order parameter A(r) = U(^f^(r)^f^(r)) and the 
molecular BEC condensate 0ju( r ) = (^( r ))- However, these are strongly coupled to each 
other and hence are not independent. One finds from Eq. (|2.1j) that their equilibrium values 
satisfy 



jj A(r) + 



2M 



0. 



This is in fact an exact identity, and can be obtained from the equation of motion 

= 



In the absence of a trap, Eq. (j2.fi j) reduces to 



9r_ 

U 



A + {2v - 2/i)0 M = 0, 



(2.6) 



(2.7) 



(2.8) 



a result discussed at length in Refs. [H, y| • As a result of this strong coupling, the Cooper- 
pair and Feshbach molecule condensates are hybridized by the Feshbach resonance, and are 
both finite throughout the superfluid phase. 



In a uniform gas, the composite order parameter appearing in Eq. (|2.4jl and defined in 
Eq. ()2.5|) can be written in the form 

A = f/ cff £< c -Pl c Pt>> (2.9) 
p 

where c pCT the annihilation operator of a Fermi atom in momentum space, and 

U eS = U + —^—. (2.10) 
2v — 2/i 

U e ff describes an effective pairing interaction in a BCS-type Hamiltonian in Eq. (j2.4j) . which 
can be tuned by adjusting the molecular threshold energy 2v. We note that the expression 
in Eq. ()2.9j) has the same form as the usual definition of a BCS Cooper-pair, apart from the 
effective pairing interaction U e g. 

The molecular Bose excitations are described by the field operator <5$(r) = $(r) — 
($(r)) = $(r) — m (r). When we substitute this expression into the boson kinetic term 
[the bilinear term involving $(r)] in Eq. (j2.1|) . the terms linear in 5$(r) are cancelled out 
by another term which is linear in <5<3>(r) arising from the last term in Eq. ()2.1j) . This is 
a consequence of the key relation in Eq. (j2.fi j) . This leaves a term in Eq. ()2.1|) which is 
bilinear in 5$(r), namely 

H M = J dr5&(r)[-^ + 2u + V* p (r)-2 f i]6*(T). (2.11) 

Equation (j2.11|) gives the lowest molecular excitation energy E^f = 2v + (3/2)co> _ 2fi. 
A self-consistent calculation shows that this threshold energy is always positive (i.e., 2v + 



(3/2)^0 > 2/x). As discussed in Ref. s| for the uniform case, this contradicts the fact 
that excitation spectrum must be gapless in a uniform interacting Bose gas (of molecules). 
Although Eq. (|2.1j) does not explicitly involve an interaction between the molecules, an 



34j. This effective repulsive 



effective repulsive interaction is induced by the Fermi gas on 
interaction will lead to a Bogoliubov phonon as the collective mode in the molecular gas in 
the BEC regime. 

In an interacting Bose gas, it is well known that the only low-energy excitations are 
collective modes. This is because the single-particle excitations are strongly hybridized |3] 
with the two-particle excitations (including collective modes). When we include the effec- 
tive molecule-molecule interaction in a consistent way, the gapless or phonon behavior of 
molecular Bose excitations in a uniform gas must be recovered. This requires an extended 



version of the present theory, as we shall discuss in a future paper j^fj. The present pa- 
per is mainly concerned with the single-particle excitations of a Fermi superfluid, and does 
not deal with the collective modes. The fact that our present treatment does not include 



the correct interaction 



33, 



34j between molecules in the BEC limit is of little importance 



when discussing the spectrum of the single-particle Fermi excitations (which disappear as 
we approach the BEC limit). 

The mean-field HFB Hamiltonian in Eq. (|2.4|) is formally identical to a trapped superfluid 
Fermi gas in the standard BCS treatment, apart from the replacement of A(r) by the (self- 
consistent) composite order parameter A(r). Bruun and co-workers have presented detailed 
numerical results for the Coo per p air order parameter and the single-particle BCS excitations 



in a BCS superfluid at T = 



37, 



38, 



jraer p£ 

3,0- 



Our present work may be viewed SIS db natural 



extension to include the effect of the Feshbach resonance, based on the model in Eq. (|2.4|) 
involving the effective pair potential A(r). The latter, of course, involves the molecular 
condensate <pM (r) and must be computed self-consistently. 

As usual, the HFB Hamiltonian in Eq. (|2.4|) can be diagonalized as -£/hfb = J2na Enlialno 
by solving the Bogoliubov de Gennes (BdG) equations |4jj, 

^ nil \ i i" i i i n .. i i" i i i ii ,. i r i i 

(2.12) 




v n r 



where H is the diagonal component of -£/hfb- The Bogoliubov quasiparticle excitations are 
described by the fermion operators j na . which arc related to the fcrmion field operator ^(r) 



as 



41] 



* T (r) 



X![«n(r)7nT + K( r hli 

n 

H[«n(r)7n| -<(r)7it 



(2.13) 



As noted above, the solutions of these BdG equations for a trapped Fermi gas have been 
discussed extensively by Bruun and co-workers 

0, y, y, q 

in the BCS limit and where 

A(r) = A(r). It is convenient to expand the fermion field operator ^(r) with respect to 



the eigenfunctions of a harmonic potential V t f (r) as 



*«r( r ) = X) fnlm( r ) C nlma- 



(2.14) 



nlm 



Here, /^ m (r) = i?^(r)Yj m (^), where Yi m {§) is a spherical harmonic and -R^(r) is the usual 



radial wavefunction, given by 



fl£(r) = V2{m^fl^ j- - t n + — e-^L^{f) (f = V^r), (2.15) 

where L^ +1//2 (f 2 ) is a Laguerre polynomial. The HFB Hamiltonian ifnFB m Eq. (|2.4jl can 
then be reduced to 

-f^HFB — E £,nl C nlma C nlma ~~ E Jnn' C nlma C n'lma 

nlm,a nn'lm,a 

~ E ^nn' \_ c nlm^ c n'l,-ml + ^- c - ■ (2.16) 
nn'lm 

Here ^ = ^o; (2n + / + 3/2) — /i are the single-particle excitation energies of the atoms in 
the harmonic potential. F l nn , and J l nn , describe the mean field effects associated with the 
composite pair potential A(r) (which plays the role of an "off-diagonal" potential) and the 
Hartree potential -|np(r), respectively. These are given by 

roc _ 

F l nn , = / r 2 drR^(r)A(r)R^(r), (2.17) 



o 



J l nn , = / r 2 drR^(r)n F (r)R^(r). (2.18) 
J o 

We note that F^ n , and J l nn , include both the intra-shell terms (n = n') as well as the inter- 
shell terms (n ^ n'). The local Cooper-pair order parameter and fermion density are given 
by0 

21 4- 1 

A(r) = Uj2^—Rni(r)R^(r)(c nm c n , m ), (2.19) 

nn'l 

21 + 1 

MO = E ^(O^(0(i^ «x>. (2.20) 

nn'lcr 

In Eq. f!2.20|) . we have taken advantage of the spherical symmetry of our model, which leads 

to (CnimiCnU-mi) = (c„,/ |<Vzot) all d ( C nlma C n'lma) = ( c n l0a C n'l0a) ■ 

It is important to remember that Eq. f)2.16j) includes the effect of the molecular conden- 
sate 0m(O' smce it enters the composite order parameter A(r) = A(r) — g T 4>M{f)- Since 
0m (0 only depends on r, we need only consider the I = quantum number. Thus we can 
expand 0Af( r ) i n terms of the radial components 

M0 = "i=E«nfl£(0- (2.21) 
V 47T n 



9 



Here f^ m (r) = R% (r)Yi m (6) is a molecular eigenfunction for the isotropic harmonic potential 



V^l p (r), with the molecular energy E^j = hu (2n + I + 3/2). The radial component R%(r) 
is identical to that given in Eq. ([2.15J1 . except that the atom mass m is now replaced by the 
bound state mass M = 2m. Substituting Eq. (j2.21|) into Eq. (|2.6|) . we obtain an expression 
for a n in terms of A(r), 



a. 



9 



UEH + 2„-2»h *^<r)i«<r). 



(2.22) 



The magnitude of the various expansion coefficients a n clearly determines to what extent 
the molecular condensate 4>m {f) in Eq. (|2.21|) is similar to the BEC order parameter 0^/ al (r) 
of a non-interacting Bose gas of molecules in a harmonic trap. The latter is given by the 
macroscopic occupation of the lowest (n = 0) state 

1 



i ideal / 
'M \ 



a Roo{r 



(2.23) 



where |ct | = yN/2 for a non- interacting Bose gas. 

Since Eq. (j2.16J) can be written as -£/hfb = Ylmi Hrfbi one ma y independently diagonalize 
each i^HTB; i- e -; f° r each set of values (I, in), by a unitary transformation. This is the 
Bogoliubov transformation, 



/ r \ 



CNilmf 
c 0l,-ml 



4/ 



iNilm] 



(2.24) 



Here W l is a 2(Ni + 1) x 2{Ni + l)-orthogonal matrix, and the fermion operator ■jnima 
describes the Bogoliubov quasi-particles. Physically, the matrix elements of W l describe the 
hybridization of particle and hole excitations in the Bogoliubov quasi-particles described by 
Inima- W l is determined by the requirement that W l ^H^ B W l be diagonal. The matrix 
elements are then obtained from the eigenfunctions of the following BdG equations for 
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/ £F — U jl 
2 J 00 



_U jl 

2 J iVzO 

-F l 



P F 



U jl 

" 2 J 07V ; 



1 J N l N l 
- F 0N t 



-Fl 



N,0 



N,N, 



-F l 



-pF 4-U T l 

SOJ ^ 2 J 00 



2 J iVi0 



-F l 



NiNi 

U jl 
2 J 0N t 



-pF + U jl I 



W l Nl+1>n 



E 



Til 



w, 



N u n 



w 



N t +l,n 



\ W 2N l+ 2,n J 



(2.25) 



As usual, we introduce a cutoff uj c = Hu!q(N c + 3/2) in the energy summation in the gap 
equation [see Eq. ()2.27|) below]. This defines N c . The maximal radial quantum number 
Ni in Eq. ()2.25|) is then given by the largest integer bounded by (N c — I)/ 2. As with the 
usual BdG equations, positive and negative eigenenergies (E^ and —E F ,) are obtained from 
Eq. (|2.25j) . The diagonalized Hamiltonian can be written using only the positive energy 
eigenenergies (E F t > 0), namely 



HFB 



E [ E nl1 ] nlm\1nlm^ ~ E^nlmllllml 



n=0 



= J2 E nl+H E F n all ma lnlm«. (2.26) 
n=0 n=0,a 

In the following discussion, we only take the positive eigenenergies E^ as in Eq. ()2.26|) . 
According to Eq. ([2.26)1 . the operator 7^ mcr describes creating a Fermi single-particle ex- 
citation from the ground state, with excitation energy E^\ > 0. Because of the assumed 
spherical symmetry of our trap, E^ and W do not depend on the quantum number m. 
Substituting Eq. (l2~2^) into Eqs. (gHHj) and QTZty . we obtain (at T = 0) 



nn 1 1 



N, 



3=0 



+l,Ni+3' 



nn'l 



21 + 1 
An 



Nl r 

'+l,Ni+j +n 

3=0 



W l - 

,j+l VV N^n'J+l 



[2.27) 
[2.28) 
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where Ni = Ni + 2 and the prime in J2' refers to the finite cutoff uo c in the summation (see 
above). 

In the BCS-BEC crossover region, we shall find that the Fermi chemical potential /i 
deviates strongly from the usual BCS limit, where it equals the Fermi energy e F . This effect 
is taken into account by considering the equation for the number of atoms in addition to the 
BdG equations^]. Since non-condensed molecules are absent at T = 0, the total number 
density of Fermi atoms n(r) is simply given by 

n(r) = 2\4> M (r)\ 2 + n F (r) = 2n u {r) + n F (r), (2.29) 

where each molecule counts for two Fermi atoms. The equation for the total number of 
atoms, which determines /i, is then obtained by integrating Eq. ()2.29|) over r. The result is 



N = 2E«n + E(^ + l)[l^ +1 ,^ +n ,| 2 + |^ 



nn'l 
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Ni+n,n'+l\ 



= 2N M + N F . (2.30) 

Here Nm and Np are the number of (Feshbach) molecular bosons and number of Fermi 
atoms in the crossover region, respectively. Equations (J2.25|) and Eq. (|2.3U|) are the basic 
equations of our theory at T = 0, taking into account a Bose condensate of Cooper-pairs 
and molecules in a self-consistent manner. We numerically solve the BdG equations ()2.25|) 
together with the generalized number equation in Eq. (|2.3U|) . determining the coefficient a n 
[see Eq. ()2.22|) ]. A(r), (fi M (r), n F {r) and /i self-consistently. 

In order to understand the essential physics, we end this Section by recalling what the 
above formalism reduces to for a uniform superfluid Fermi gas. In this case, the simple pair- 
ing mean-field approximation (MFA) gives usual BCS-Gor'kov expressions for the diagonal 
and off-diagonal single-particle Green's functions 

G n (p,u) = ^ 2 + _% ' Gl2(p,w) = "^T^' (2 - 31) 
with the BCS-Bogoliubov excitation energy given by 



E P = V(£ P -^) 2 + |A| 2 . (2.32) 

Here £ p = e v — fi (where e p = p 2 /2m) is the kinetic energy measured from the chemical 
potential. Using G±2(p,uj) to calculate A, one finds that A in Eq. ()2.9j) satisfies the usual 
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"gap equation" [but now with the pairing interaction U e g given by Eq. (|2.10J) ]. 

A = U ^Y,'^r tanh \(3E P , (2.33) 

valid at finite temperatures. As usual, a cutoff is introduced in the momentum summation. 
At this MFA level, the total number of atoms at T = has the same form as Eq. (|2.30p . 

N = 2N M + N F . (2.34) 

The number of Fermi atoms in Eq. (|2.34|) is now given [using Gn(p,w)] by the well-known 
BCS expression 

= £<4, CT c pa > = El 1 - It tanh \f3E± (2.35) 

while the number of condensed molecules is Nm = \4>m\ 2 - 

In the uniform gas, the values of \i and A are determined by the self-consistent solutions 
of the MFA gap equation ()2.33|) and the number equation given by Eqs. ()2.34|) . This 
simple "pairing approximation" for the single-particle Fermi excitations is expected to give 
a quantitative description at T = [where tanh \(3E P — > 1 in Eqs. ()2.33j) and (|2.35j) ]. since 
fluctuations are small and all the molecules are Bose-condensed. This T = limit was first 
studied by Eagles[8] and Leggett[9] in the absence of a Feshbach resonance. As we mentioned 
earlier, for T approaching T c , the fluctuations associated with exciting molecules out of the 
condensate and coupling to the particle-particle (Cooper-pair) channel become dominant. 

These fluctuations (rather than the breaking up of two-particle bound states) were first 
included by Nozieres and Schmitt-Rink (NSR)j3| to determine T c . In Ref. jsj], we extended 
this NSR approach to discuss the superfluid phase below T c in a uniform Fermi gas, including 
a Feshbach resonance and associated molecular bosons. Both A and /i are obtained by 
solving Eqs. ()2.33|) and ()2.34|) self-consistently, but now include the depletion of A through 
the presence of non-Bose-condensed molecules. This procedure gives the simplest extension 
of the MFA-BCS single-particle results, in that now Eqs. (j2.31)l and ()2.32|) involve values 
of \x and A which include (in an average way) the effect of fluctuations around the MFA 
theory. The effect of such fluctuations is of considerable current intere st 143 . 
what is called the "pseudogap" regime (for a recent review, see Ref. |45|). In this region, 
strong low-energy fluctuations in the Cooper-channel suppress the density of states around 
the Fermi energy, which has the same effect as if there was an effective pair potential, even 



44j and lead to 
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outside the superfluid region (T > T c ) where A(r) vanishes. In a future paper, we will use 
this generalized NSR approach to include fluctuations and extend the results of the present 
paper to finite temperatures. 



III. SINGLE-PARTICLE GREEN'S FUNCTIONS 



In ordinary BCS superfluidity, the single-particle excitations (BCS quasiparticles) are 
associated with dissociation of weakly-bound Cooper-pairs. In the coupled fermion-boson 
model in Eq. (|2.4jl . these Cooper pairs are replaced by composite bosons, consisting of 
Cooper-pairs and molecular bosons associated with the Feshbach resonance \A(r) = A(r) — 
<?r0M( r )]- Even in the BEC regime (2u < 0), Fermi excitations can still exist [30] as well 
as collective modes (which form a Bose spectrum). As more and more fermions pair up to 
form bosons, the spectral weight of the Fermi branch vanishes, shifting to the Bose collective 
branch. 

Single-particle properties are most conveniently discussed in terms of Green's functions. 
The "diagonal" single-particle thermal Green's function Gn(r, r', iu m ) is defined by 



48] 



47 



G n (r,r',iu; m ) =-jyre ib *» T (T r {*i(T,T)*\(T>,0)}), (3.1) 

where iu m is the fermion Matsubara frequency associated with the imaginary time r. One 
needs three other single-particle Green's functions to describe a Fermi superfluid, as sum- 



7] 



marized by the 2x2 matrix Green's function [4 

Gu(v,r',iuj m ) G 12 (r,r',iuj 
G 2 i(r,r',iuj m ) G 2 2(r, r',iu. 



G(r, r', iu % 



in J 



- / dre luJmT 



o 



T T {* T (r,T)*{(r',0)}}, (T T {* T (r,r)^(r',0)} 
T r {*!(r,T)*{(r',0)}}, (T r {*I(r,r)^(r',0)} 



(3.2) 



In Eq. ()3.2)1 . G22 gives the single-particle excitation spectrum of Fermi atoms of pseudospin 
|. The off-diagonal components Gu and G21 arise as a direct consequence of the broken 
symmetry and the presence of a condensate of Cooper pairs. Using the BdG equations 
in Eq. ()2.12j) . one can show that these Green's functions are related to each other as 
G r 2 2(r / , r, ^m) = — Cn(r, r', iu m ), and G^r', r, — iu m ) = Gi 2 (r, r', iu m ). From the definition 
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of the Cooper-pair order parameter A(r) = f/(\&|(r)\I / |(r)), we find the important self- 
consistency condition (the gap equation) 



A(r) = ±£'G 12 (r,r, 







(3.3) 



Using the eigenfunctions /^ tt (r) defined in Eq. (|2.14|) . we can write 2x2 single-particle 
Green's function in Eq. (|3.2|) as 



G(r,r',iu m ) =J2Yi m 0)g l (r,r',iu m )Y l * rn 0'), 

lm 



(3.4) 



where g l (r, r', iu m ) is the 2x2 Green's function for a given value of the angular momentum 



g l (r, r',iu m ) 



E 

3=0 



A i ,(r)4(r') , A,,(r)AWr' 



+ 



iuj m — En iuj m + En 



(3.5) 



The two-component spinor Aji(r) in Eq. ()3.5|) is defined in terms of the solutions of the 
BdG equations in Eq. ()2.25|) . namely 



Aii(r) = E 



N \ ( W l 



n=0 



w l 

" + J -' +J 1 'Or)- A ;/ (r) *Y 2 A,,(r). 



(3.6) 



Here r 2 is the Pauli matrix. 

A very useful quantity is the Fermi single-particle excitation spectrum N(u), also referred 
to as the density-of-states (DOS). This is related to the spectrum of single-particle Green's 
function, 



7T 



N(u>) = / dvlm Gn(r, r, iu m — > u + iS) 



(3.7) 



Substituting Eqs. ()3.4j) and ()3.5|) into Eq. (|3.7|) . this spatially-averaged density-of-states 
(per pseudospin) is given by (En > 0) 



N(u) = E(2/ + 1) E \K +W \ 2 ^ - Eft) (u > 0). 

I n,j 

It is also useful to introduce the local density of states (LDOS), defined as 

N(u,r) = Im Gn(r, r, icu m — > u + iS) 

91 4- 1 

- E -^Rni{r)Kn{r) E W l n+w W l n , +w 5(^ - Ef>) 



(3.8) 



nn' l 



(3.9) 
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The local density of states LDOS in Eq. ()3.9|) is simply related to the spatially-averaged 
total DOS in Eq. by 

N(u)= A7rr 2 drN(uj } r). (3.10) 
Jo 

The gap equation (|2.27|) is easily obtained from the present Green's function formalism. 
When we use the off-diagonal (l,2)-component of Eq. ()3.5|) in Eq. (j3.3|) . we obtain 

A(r) = -UY!^-Ri {r)R F nll (r) £ W^ +nJ+1 <' + i,; + i- (3.11) 

nn'l 47F j=0 

Once one has A(r), one can calculate the coefficients a n in Eq. ()2.22|) and then find </>m(?")- 
By this procedure, one finally obtains the value of A(r), with the total number equation 
in Eq. (J2.3U)) determining the self-consistent values. As with the usual BdG equations in 
a uniform BCS model, there is a relation between eigenfunctions for E^i > and < 0, 



which is given by 



41] 



E 



w 



( 



Ni+n,j+l 
"^n+l,j+l 



Rni(r) 



E 



w. 



n+l,Ni+j 

wL 

\ vv Ni+n,Ni+j 



KM. 



(3.12) 



Substituting Eq. ()3.12|) into Eq. ()3.11|) . we find that the gap equation ()2.27j) for the Cooper- 
pair order parameter A(r) is reproduced. 



IV. SCATTERING LENGTHS, GAP EQUATION AND STRENGTH OF FESH- 
BACH RESONANCES 

The Hamiltonian in Eq. ()2.1)1 involves the bare energies U, g r and 2v. As usual in dealing 
with ultracold atomic gases, it is convenient to work in terms of renormalized interaction 
energies which incorporate the effect of high energy processes. This procedure naturally 
leads to the two-body scattering length a 2 ^ which describes the effective interaction between 
low energy atoms, even in the case of a Feshbach resonance. The two-body scattering length 
a 2 . b can be measured directly in a variety of ways. 

In this section, we briefly review the standard theory for renormalized low energy param- 
eterS (see, for exa^nple, See. IV A of Ref. D and Rets. QQQ). We also point ont that 
in the presence of a Feshbach resonance, the self-consistent gap equation which determines 
the order parameter naturally introduce a different s-wave scattering length a s , which has 
a crucial dependence on the Fermi chemical potential fi. As a result, in dealing with the 
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BCS-BEC crossover in Fermi superfluids, it seems most natural to treat a s as the control 
parameter, rather than af. 

There is a second reason which makes it useful to discuss properties in crossover region as 
a function of a s . Before doing so, we discuss the parameters we use in this paper. We take 
the total number of atoms to be N = 10, 912 (= N^ + N^ = 2N^). In a non-interacting Fermi 
gas, this corresponds to filling atoms (per spin) up to E = 31.5huo (= £f), where uq is the 



trap frequency. As the unit of length, we use the Thomas- Fermi radius R F = y2ep /muiQ 
for a free Fermi gas in the trap. We take the Feshbach coupling constant g r = 0.06u 



(g T = g I (yN/Rp) = 0.2£f). The non-resonant pairing interaction values we use are either 
U = O.OOlcuo (U = U{N/Rl) = Q.35e F ) or 0.0015cj (U = 0.52e F ). For the high-energy 
cutoff, we take u> c = 161.5^0 (3> £f)- 

The explicit calculations presented in this paper are for what is called a narrow Feshbach 
resonance, while all current ultracold Fermi gas experiments are done using broad Feshbach 
resonances. However, in a uniform Fermi gas, several quantities are found to have values very 
similar if viewed as a function of a s , for both weak (g r < e F ) and strong (g T > £p) Feshbach 
resonances. At the present time, numerical calculations are only able to deal with a narrow 
Feshbach resonance in the case of a trapped Fermi gas. This is due to the fact that one is 
limited to dealing with a finite number of excited states in a trap and this makes it difficult 
to deal with a broad Feshbach resonance, which couples molecules to Fermi atoms in very 
high-energy eigenstates. This approximate independence of the strength of the Feshbach 
resonance (for a given value of a s ) is thus very useful. 

We first recall the standard case of two Fermi atoms interacting in a vacuum[48] with a 
bare interaction denoted by —U. In this case, the effective low-energy (u = 0) renormalized 
s-wave scattering length is given by 

- ^£ = U - r = U*. (4.1) 

Effectively the renormalized interaction U R incorporates all high energy scattering processes. 
In the absence of a Feshbach resonance (g T = 0) and in a uniform gas, the gap equation 
(I2.HHJ) can be written in terms of this renormalized interaction as follows: 

l = U RJ2\J_ tanh ^ - — 1 . (4.2) 
^ L2E 2T 2f -I 

Here we have written the gap equation at finite temperatures. The summation in Eq. (|4.2|) 
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converges, so we can safely takes u> c — > oo. This gives a cutoff-independent gap equation as a 
function of the renormalized interaction U R or, equivalently, the two-body scattering length 
af defined in Eq. (|4.1|) . As first discussed by LeggettQ and Randeriajl2|. it is convenient 
to discuss the BCS-BEC crossover region in terms of the dimensionless parameter (kpaf')^ 1 . 
As the bare attractive interaction U increases, (/cpc^ 6 ) -1 goes from — oo (BCS) to oo (BEC). 

In the case of a Feshbach resonance, one can introduce a renormalized s-wave scattering 
length describing low energy atoms which is the analogue of Eq. (j4.1|) . namely, 



4rf U. 



21 j 



— | i -26 1 ' (4-3) 

m 1 - t7 eff 2^[ ,a; c ] 2iT 



where the bare Feshbach resonance is described by 



V% = U + £ (4.4) 

Once can rewrite Eq. (|4.H|) in the form 

4tTO 26 (o R ) 2 

_^_ = U R + ^L (4 .5) 
m 2u u 

which now involves the renormalized parameters for U, g r , and 2v (for details, see, for 
example, Sec. IV A of Ref. Q]). These low energy renormalized parameters can be directly 
measured. 

However, in contrast to Eq. (|4.2j) . the cutoff- independent energy gap equation in the 
presence of a Feshbach resonance in not given simply in terms of the usual two-body s-wave 
scattering length as defined in Eqs. (j4.H|) and (|4.5jl . Instead, one finds jl|, Q, 0, 0| 

l = Ujk Y [J-tanh^-— 1, (4.6) 
cff rn ^ \2EL 2T 2e p J' v ; 

where 

T R _ ^eff 



U. 



1 - Ucs E[0,u, c ] 2^ 



2u R - 2/2 

^ (4-7) 



m 

U e ff is defined in Eq. ()2.10|) . We first restrict the following discussion to the case of a 
uniform Fermi superfluid. Comparing this new s-wave scattering length a s involved in the 
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gap equation ()4.6|) with the the two-body scattering length of defined by Eqs. ()4.3|) and 
()4.5|) . we see that a s depends crucially on the Fermi chemical potential /i. This in turn is 
a strong function of the threshold energy 2v and hence will change value in the BCS-BEC 
crossover. One may think of a s as including many-body effects related to the coupling of 
the fermions to molecular bosons 3] . We note in passing that it is wrong to simply use af 
[given in Eq. ()4.5|) ] in the gap equation ()4.6|) in the general case. 

The fact that the cutoff-free gap equation ()4.6|) only depends on a s , irrespective of detailed 
values of the non-resonant interaction U and the Feshbach resonance strength g r , is very 
useful. In particular, a broad Feshbach resonance (g T > Ef) and a narrow Feshbach resonance 
(g r < e r ) give approximately the same values for T c and A for the same value of a s . As 
an example, we show in Fig. ^a) the phase transition temperature T c in the BCS-BEC 
crossover for a uniform gas for a broad and a narrow Feshbach resonances. In a broad 
Feshbach resonance, the BCS-BEC crossover occurs at the value of v ^> e-p. In the case of 
a narrow Feshbach resonance, the crossover occurs at v ~ e^. However, when we express 
T c as a function of (kpag)' 1 , both cases give almost the same result over the whole BCS- 
BEC crossover region, as shown in Fig. ^a)[50]. In addition, Fig. E^b) shows that the 
change from a gas of Fermi atoms to a molecular Bose gas in the BCS-BEC crossover is 
also almost the same for both large and small values of g T when viewed as a function of 
(kpOs)^ 1 . The region of (kpOs)^ 1 > can thus be regarded as a BEC of molecular Bosons 
for both broad and narrow Feshbach resonances. Using a s in Eq. ()4.7j) . we find the same 
crossover physics irrespective of the width of the Feshbach resonance. This is useful since as 
we noted earlier, our numerical calculations are limited to g r < e-p, while recent experiments 
deal with a broad Feshbach resonance. We also note that in the case of a broad Feshbach 
resonance, the crossover occurs at a value of v ^> ep > fi. In this region, we can omit 
the chemical potential in Eq. ()4.7|) and then one finds a s ~ af , the usual renormalized 
low-energy two-body scattering length which can measured directly. 

In contrast to T c , the character of the Feshbach resonance (broad or narrow resonance) 
does show up somewhat when we consider the number of molecules Nm in Fig. ^b). With in- 
creasing (/cpa s ) _1 from the BCS side, Fig. I2fb) shows that the number of Feshbach molecules 
with a finite lifetime increases in the case of a narrow Feshbach resonance. Stable molecules 
(Nft~°) become dominant in the region of negative chemical potential, /z < (see Fig. Efa)). 
Although stable Cooper-pairs given by the iV^ -0 curve also exist in this region, we see that 
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iV2 =0 < N^ =0 . Thus, in the case of a narrow Feshbach resonance (small g r ), superfluidity 
in the region of fi < is largely associated with a BEC of stable Feshbach molecules. 

On the other hand, in the case of a broad Feshbach resonance, the number of Cooper- 
pairs with a finite lifetime (Nq) first rapidly increases as one increases (/cpa s ) _1 [see Fig. 
EJc)]. When /i < 0, stable long lived Cooper-pairs (Nq~°) become the dominant bound 
states. Below T c , the tightly-bound molecules around /i < are now Cooper-pairs. As 
shown in Fig. Efc), the Feshbach molecules (iV^ -0 ) eventually become dominant deep inside 
the BEC regime (defined by (kpOs)" 1 ~ 1). The difference between the two cases is due the 
fact that, in a broad Feshbach resonance, the BCS-BEC crossover occurs around v ^> e-p, 
where the formation/dissociation of Feshbach molecules only appears as a virtual process. 
In contrast, the crossover region is located in the region v < ep in the case of a narrow 
Feshbach resonance, where a large number of Feshbach molecules can form as quasi-stable 
entities. 

The preceding discussion was limited to a uniform Fermi superfluid. A discussion on T c 
and the composite order parameter in terms of the scattering length a s is also possible in a 
trapped gas, if we use the LDa|2.1o"i1|. In this case, the spatially-dependent order parameter 
A(r) is self-consistently determined by the cutoff-dependent gap equation, given by 

1 = ^w [ £ ] 4w tanh ^ 1 ' (4 ' 8) 

Here E p (r) is the energy of a BCS excitation with density and other quantities evaluated 
locally at position r. This is given by Eq. ()2.32|) . with fi being replaced by fi(r) = fi — 
muo^r 2 /2. The bare pairing interaction strength at r in Eq. (j4.8|) is defined as 

U cS (r) = U+ f -— , (4.9) 

2u + (3/2)u;o — 2/i(r) 

where we have included the zero point energy (3/2)^0 for later discussions. From our 
preceding analysis, the LDA gap equation ()4.8|) can be written in a cutoff-independent form, 
but it now involves the renormalized pairing interaction given by 

U&r) = ^ -. (4.10) 

In particular, since /i(r = 0) = fi, U^{r = 0) equals given in Eq (|4.7|) [apart from 
the zero point energy term (3/2)cj i n Eq- (|4.9j) ]. Thus, recalling that T c in the LDA 3| is 



20 



determined by the gap equation (|4.8|) at r = 0, we conclude that the crossover behavior of 
T c and A(r = 0) at the center of the trap can be described (as in a uniform Fermi gas) as 
functions of the renormalized scattering length a s defined in Eq. (|4.7|) . Indeed, Fig. Efa) 
shows that both narrow and broad Feshbach resonance cases give almost the same T c as a 
function of (/cf^s) -1 in the whole BCS-BEC crossover regime. 

As in a uniform gas discussed above, in a trap some quantities will depend on the strength 
of the Feshbach resonance. The LDA gap equation ()4.8|) involves a position-dependent a s (r) 
because of the spatial dependence of //(r). In a broad Feshbach resonance (g T ^> Ep), the 
crossover occurs at v sy, so that /i(r) (~ ep) can be neglected in this regime. On the 
other hand, U e s(r) decreases as we go from the center of the trap in the case of a narrow 
Feshbach resonance (g r < £p), since now the crossover occurs at v ~ e-p and hence /x(r) 
cannot be ignored in U e ^{r). This explains why the total number of molecules in slightly 
smaller in the case of a narrow Feshbach resonance, as shown in Fig. Efb) 

In this paper, we present results as a functions of the renormalized uniform gas scattering 
length a s in Eq. (J4.7j) . with the result that the theory is approximately cutoff-free within 
LDA. The dependence on (/cp^s) 1 also gives a description for Fermi superfluids around 
the center of the trap, irrespective of the width of the Feshbach resonance. We refer to the 
region of (kpCis)^ 1 < as the BCS regime and the region of (kpag)^ 1 > as the BEC regime. 
We note that the "real" or bare pairing interaction in the gap equation is not 47ra s /m, but 
— U e g defined in Eq. 1)2.10)1 . This is always attractive even in the positive (BEC) region of 
a s > 0, as we show in Fig. Ufa). In Fig. IHb), we give the relation between a s and z/, for the 
parameters g r and U used in this paper. In the numerical results discussed in this paper, 
for the convenience of the reader, we give both the value of v as well as the corresponding 
value of (kpas)^ 1 . 

V. EQUILIBRIUM PROPERTIES IN THE CROSSOVER REGION 

Figure El shows the calculated atomic density profile at different places in the BCS-BEC 
crossover region at T = 0. In the BCS regime shown in Fig. Efa), the number of Fermi 
atoms is much larger than the number of (Feshbach) molecules. Although these atoms form 
Cooper-pair bosons, the Pauli exclusion principle is still relevant. As a result, the Fermi gas 
spreads out to a distance of the order of the Thomas- Fermi radius Rf, just as in the case of 
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a trapped non-interacting Fermi gas. We recall that in contrast to bosons, interactions have 
a negligible effect on the density profile of a trapped Fermi gas in the normal phase. The 
relative number of Cooper-paired fermions [wp(r)] and bound (Feshbach) dimers [nn(r)] 
continuously changes as the threshold energy 2v is lowered (or (/cpa s ) _1 increases). On 
the BCS side of the crossover regime shown in Fig. Efb), the molecular density %(r) has 
become very large in the center of the trap. This shows that the Pauli exclusion principle, as 
expected, is less important in this region where Bose molecules start to dominate. However, 
the Fermi atoms still dominate for r > 0.2Rp. This feature persists on the BEC side of 
crossover regime shown in Fig. E3[c), where np(r) is seen to be dominant for r > 0.3i?F- The 
tail of the density profile consisting of unpaired atoms finally disappears in Fig. E(d), where 
almost all the atoms have formed Feshbach molecules. 

Figure El shows how the Fermi chemical potential fi changes in the BCS-BEC crossover. 
As one approaches the BEC regime, \x decreases from the weak-coupling (BCS) result given 
by \i ~ ep (= 31.5co>o in our case). In the BEC regime, the molecular Bose chemical potential 
Hm = 2/i approaches the lowest molecular energy, given by 2v + (3/2)u;o, as shown in the 
inset in Fig. El We note that this limiting case of \im = 2u + (3/2)cl> is just the condition 
for the BEC of a non-interacting Bose gas in a harmonic trap. 

In Fig. the Fermi chemical potential is seen to go through zero at (kpa s ) _1 ~ 0.65, 
which corresponds to v — 0. To understand this, let us consider the case of a uniform gas. 
In this case, the single-particle excitation spectrum is given by the usual BCS expression 
(but now with the composite order parameter A) P. I12I l^. 




In the weak-coupling BCS regime, where fi is positive, the energy gap of the BCS excitations 
at the Fermi energy is equal to |A|. However, this relation is no longer valid when /1 is 
negative, as in the strong-coupling crossover regime 0, Q . In this regime, the energy gap is 
given by 



In the BEC limit {y — > —00), the composite order parameter is described by a pure BEC 



condensate of molecules, given by A — >• g r ^jN/2. In the same limit, the chemical potential 
approaches the threshold energy \i ~ v — ► — 00. Thus, the excitation gap E g given by Eq. 







(5.2) 
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()5.2|) approaches \u\. This makes sense, since 2u (< 0) is the lowest molecular energy, and 
hence 2\v\ is the excitation energy to dissociate a molecule into two Fermi atoms. 

In a trapped gas, within the LDA, the chemical potential // and the order parameter A 
in Eq. (|5.2|) are replaced by the position-dependent ones, fi(r) = fi — muj^r 2 /2 and A(r), 
respectively. In the extreme BEC limit (v — > — oo), the lowest single-particle excitation 
energy is given as the energy to dissociate a molecule into two Fermi atoms and put them at 
the lowest (unpaired) fermion state. In the LDA, this lowest state is at r = 0, because the 
trap potential has its minimum at the center [VjJ (r = 0) = 0]. This again leads to E g ~ 

~ |z/| in the BEC limit (when A(0) <C | /x | ) . As a result, we find that {k-pa s )~ l ~ 0.65 (at 
which fi = 0) gives a characteristic scattering length which separates the region dominated 
by weakly-bound Cooper-pairs (fi > 0) and the region dominated by tightly-bound molecules 
(/i < 0). Equivalently, this boundary occurs at v — 0. 

Figure [7| shows the profile of composite order parameter A(r). As expected from Fig. 
A(r) is dominated by the Cooper-pair component A(r) in the BCS regime shown in panel 
Ufa). This calculated profile agrees with previous results in the BCS limit, as obtained 
by Bruun and co-workers [^3]. However, in the presence of a Feshbach resonance, even in 
the BCS regime, the molecular condensate <pM{ r ) is found to be finite everywhere A(r) is 
finite, due to the identity given in Eq. ()2.6|) . In an ideal Bose gas BEC, </>iw(r) is simply 
proportional to the ground state wavefunction of the harmonic potential given by Eq. (j2.23|) . 
However, since the molecular bosons are strongly hybridized with the Fermi atoms by the 
Feshbach resonance, in the BCS regime, we find that 0a/ {r) is no longer given by the ground 
state of the harmonic potential as in Eq. ()2.23j) . Indeed, as shown in Fig. IHfa), 0m(^) in 
this regime is given by a superposition of excited molecular states (n > 1), as indicated in 
Eq. (j2.21|) . The molecular condensate 4>m{^) continuously increases in magnitude as one 
approaches the BEC regime, as shown in Figs. E^b) and (c). At the same time, Fig. |H] shows 
that the contribution of excited states (n > 1) decreases, indicating that 4>m(j) approaches 
the ideal BEC described by Eq. fl2.23j) . However, because of the identity in Eq. (|2.fi|) . the 
Cooper-pair component A(r) remains finite even in the BEC regime. Indeed, as shown in 
Fig. [7(d), A(r = 0) ~ IbujQ when v = 0. 

In Fig. El we compare the total atomic density profile n(r) = 2n^i(r) + n-p(r) with the 
profile of the composite order parameter A(r), which clearly shows that the composite order 
parameter is always finite in the region where n(r) is finite. In addition, we expect that 
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n(r) ~ \4>M{f)\ 2 in the BEC regime. This is confirmed by the calculated values shown in 
the inset in Fig. E^b), which show that n(r) ~ |A(r)| 2 . 

We note that when the LDA is used in calculating the density profile as well as the order 
parameter at T = 0, it has been shown that the LDA is a good approximation in the BCS 
regime In contrast, the shrinkage of the density profile n(r) shown in Fig. as well as 
the profile of the composite order parameter A(r), is poorly overestimated by the LDA 51 1 
as we enter the BEC regime. This is because the LDA underestimates the kinetic energy, 
which results in the density profile of atoms spreading out more. 

At T = 0, such an LDA calculation correctly predicts that the resulting composite order 
parameter is finite in all regions where n(r) is finite. This reflects the fact that the Fermi 
surface is unstable against an infinitesimally weak attractive interaction at T = 0, which 
leads to a superfluid phase transition everywhere (as long as n(r) is positive) in a trap in 
the LDA. However, this situation is no longer satisfied at finite temperatures. As we have 
discussed in Ref. |2j , only the center of the trap (r = 0) is in the superfluid phase just below T c 
in the LDA, with this superfluid region becoming wider with decreasing temperatures. Thus, 
for < T < T c , an LDA calculation predicts a spatial region where the order parameter 
vanishes, even though the particle density is still finite. (This feature is shown in Fig. 1 of 
Ref. |29j.) However, the presence of such a "two-phase" trapped Fermi gas (superfluid at 
the center and normal phase at the edge of the trap) is clearly an artifact of the LDA. In 
fact, the entire gas is in the superfluid phase below T c , with A(r) finite everywhere where 
n(r) is finite. 

VI. SINGLE-PARTICLE EXCITATIONS IN THE CROSSOVER REGION 

Figure ITU1 shows single-particle density-of-states N(u) given in Eq. (|3.8|) . One sees that 
the excitation spectrum in a Fermi superfluid has a finite energy gap (= E g ). Although the 
pairing interaction becomes stronger with decreasing threshold energy 2z/, the magnitude of 
E g ~ 1.2cg> is almost the same in panels (a) and (b) of Fig. We also note that there is 
characteristic sharp increase or shoulder around E g . Recalling the case of a uniform BCS 
superfluid (where the energy gap is equal to |A|), one has the well-known result 



N(u) = p(e F ) 




e(cu- |A|) 



(6.1) 
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where Q(x) is the step function. The normal-state density of states is given by p(u) = 
. [We have approximated this by the value at the Fermi energy, uj = ep, in Eq. (j6.1|) . 
This is a good approximation in the BCS limit. The more correct expression for N(uj), 
including the effect of the energy- dependence of p(uj), is given in Ref. [21].] Equation (j6.1j) 
predicts that N(u) is singular at to = |A|. The sudden increase in N(u) at uj = E g in the 
trapped gas shown in Figs. ITUT a) and (b) may be viewed as the remanent of this singular 
behavior of the single-particle excitations in the uniform BCS case. 

Figure ITUT c) shows the change in the density of states N(uj) as we enter into the BEC 
region {y < 0). One finds that the energy gap now sharply increases (E g ~ 5u>o) and the 
sudden jump in the magnitude of N(uj) at the energy gap E g is absent. 

7 ig . ITT] we plot the single-particle excitation gap E g in the BCS-BEC crossover 
53[. This figure shows that E g rapidly increases when the chemical potential p be- 
comes negative [(^f^s) -1 > 0.65]. In this region, 2E g is the dissociation energy of a molecule. 
This energy 2E g ~ \pm\ — 2\u\ becomes large as we decrease the threshold energy 2v — > — oo. 
Indeed, we find that E g approaches \p\ (given by the dashed line in Fig. [TT]) when we take 
p — > — oo, which is consistent with Eq. (15. 2j) valid for a uniform gas. We also see from Fig. 
IT21 that the number of molecules becomes much larger than the number of Fermi atoms for 
v < or (/c F a s ) _1 > 0.65, as expected. 

One can understand (in the case of a uniform gas) why the "coherence peak" (i.e., the 
sudden jump) in N(uj) at E g is absent in the crossover region shown in Fig. fToT c). In the 
BEC region one finds (see Eq. (5.4) of Ref. j^j]) 

N(u) = ^^(l + U , ) ( y/u* - A 2 + p) l/2 e(u - (6.2) 



In 

region 



where the factor ( y u; 2 — A 2 +p) 1 / 2 comes from the normal state density-of-states p(u) oc y/uj. 
When p is negative, the threshold energy of the N(u) is determined not by the factor 



to 2 — A 2 , but by the threshold energy of the normal DOS p(uj). This leads to the step 
function in Eq. ()6.2|) . As a result, the density-of-states N(u) of a Fermi superfluid is finite 



only for to > yp 2 + A 2 . The expected coherence peak at uj = |A| < y p? + A 2 actually 
occurs in a region where N(uj) vanishes. 

When we compare the single-particle excitation gap E g for p > in Fig. ^2 with the 
magnitude A(r = 0) of the composite order parameter shown in Fig. El we see that E g 
|A(0)| ~ ep in the crossover region. Since the excitation gap E g equals the order parameter 

25 



(when fi > 0) in a uniform gas, this difference clearly originates from the inhomogeneity of a 
superfluid Fermi gas in a trap. To understand the origin of this very small value of the single- 
particle excitation gap E g in a trapped superfluid Fermi gas, it is very useful to compare N(cu) 
and the local density of states N(u,r) given by the expression in Eq. (J3.9|) . Figure ITlT b) 
shows that no low energy spectral weight comes from the center of the trap. In fact, panel 
(d) clearly shows that the low-energy spectral weight in N(u) shown in panel (a) mainly 
comes from the region r ~ 0.6-Rf, which is close to the bottom of the effective potential well 
composed of the (composite) pair potential A(r) and the harmonic trap measured from the 
chemical potential /i, given by V^L(r) = [V t f ap (r) — A0@(Kra P ( r ) — jf) [see the inset of Fig. 
14(b)]. Strictly speaking, this effective potential is a combination of an ordinary (diagonal) 
and an anomalous (off-diagonal) potential in the BdG equations. This situation is analogous 
to the boundary problem in superconductivity, as schematically indicated in Fig. As first 
discussed by de Gennes and Saint- James (o^ in the context of superconductivity, "Andreev" 
bound states can appear well below the bulk energy gap (~ A), around the minimum of 
the effective potential well. Thus E g may be viewed as the energy of the lowest Andreev 
bound state formed in the effective potential well [= A(r) + V^L(r)]. These are also called 



trap 

the "in-gap" states[57|. 

The role of such surface excitations was first discussed by Baranov|57j in connection 
with a BCS superfluid gas in a trap, using the WKB semi-classical solution of the BdG 
equations. More recent papers in the context of the BCS-BEC crossover are by Kinnunen, 
Rodriguez, and Torma using LDA, and by Heiselbergj^] using the WKB approximation. 
We note that bound state energies decrease when the effective potential width d in Fig. ITHTb) 
increases. In a trapped Fermi gas, the analogous effective potential well becomes wider as we 
enter the crossover regime because the spatial width of the composite order parameter A(r) 
shrinks in this region, as shown in Fig. [7| In addition, the decrease of chemical potential 
\i leads to a more gradual slope of the diagonal potential V^ p (r) around V^f p (r) = 0. 
These are the reasons why E g slightly decreases (see the inset of Fig. ITT]) with increasing 
(kpag) -1 , in the region 0.65 > (kpag)^ 1 > 0. We note that the lowest Andreev bound state 
energy level does not determine E g once we are in the BEC region (fi < 0). In this case, 
2E g is dominated by the dissociation energy of tightly bound molecules and thus we find 
E ~ I M I ~ M (v — > — oo) in Fig. [TT] 

nn 

In connection with the WKB approximation|57l. 1581] . we remark that this should be quite 
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good in the extreme BCS limit, where there is still a well-defined Fermi surface. However, 
this approach seems less justified in the crossover region, where the important single-particle 
states are no longer close to the Fermi energy. Further studies are needed of the validity of 
the semi-classical approximation to the BdG equations in the crossover region. 



VII. RF-TUNNELING SPECTROSCOPY IN A TRAPPED FERMI SUPER- 
FLUID 



In this section, we discuss recent work using rf-tunneling spectroscopy. As with any 
tunneling experiment, this clearly gives information about the single-particle excitations in 



fermion superfluids[27 



29 



54, 



55 



56 



58, 



59 



60(. In particular, one can extract information 



about both the energy gap of the single particle excitations as well as the value of the pair 
order parameter, but this can only be done by comparison with theoretical calculations. 
The rf-tunneling spectroscopy cross-section is calculated by considering the tunneling current 
induced by laser radiation. In the rotational wave approximation, the tunneling Hamiltonian 



is given by [29, 



54 



55 



56] 



= t F f dr[e i(qL - r " WL ' } *i(r)^ T (r) + h.c. 



+ t 



M 



dr 



(7.1) 



Here and ujl represent the momentum and frequency of the laser light, respectively. The 
first term Hf describes the usual tunneling of an atom in the j-spin state to another hyper- 
states described by the fermion field operator ty a (r), with a strength given by the tunneling 
matrix element tp. (We assume that tp has no spatial dependence.) This new hyperfine 
state \l/ a is described by the Hamiltonian 



H„ 



V 2 

-2^ + ^-Ha + V t 



F 

trap 



(7.2) 



where u a is the threshold energy and fi a is the chemical potential of this state. In Eq. (|7.2j) . 
we assume that atoms in this state are non-interacting and feel the same trap frequency as 
the other hyperfine states |, j. The second term Hf 1 in Eq. (|7.1|) describes the tunneling 
into state \l/ a associated with the dissociation of a bosonic molecule, with the matrix element 
t M . This process is the signature of a Feshbach resonance. In the superfluid phase at T = 0, 
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the Bose quantum field operator $(r) can be replaced by the macroscopic wavefunction 
0jw(r). In this case, the molecule component in Eq. (j7.1|) reduces to 



H t M = t M [ dr[e l ^™'V A f(r)^(r)^J(r) + h.c 



(7.3) 



The tunneling current operator is obtained from I(t) = N a (t) = i[H, N a (t)], where N a 
j drfyj(r) ty n { r) is the number operator of the \l/ a -state (for details, see, for example, Chap. 
9 of Ref. |46j). The resulting current operators originating from Hf and H^ 1 are given by, 
respectively, 



I F (t) = -it F / drle^^-^^Kr^^r) - h.c. 



i M {t) 



-it 



M 



dr 



e i(qL -™'VAf(r)^i(r)^I(r) - h.c. 



(7.4) 
(7.5) 



Assuming the tunneling matrix elements t F and tM are small, we can evaluate the tunneling 
current using first order perturbation in H t . The current associated with the usual tunneling 
term Hf induced by the rf-field is given by (we set = 0) 



I f (lu) = (I f (uj)} = 2t F Im / drdr'H F (r, r', — uo 



where 



U) = LO L - LU a - fl + fl a 



(7.6) 



(7.7) 



defines the effective detuning frequency. The two-particle Green's function n^(r, r', u) in 
Eq. ()7.6)1 is obtained from the analytic continuation)^] of the thermal Green's function 



U F (r,r',iiy n ) 



dTe iv » T (T T {*l (r, r)¥ T (r, t)*Ut>)V a (r')}) 



~a G ll( r ' r '' iuJ m + iv n )G a {r', T, lU m ) . 

P iaj™ 



(7.8) 



Here, G a is the single-particle thermal Green's function for the \l/ a state quantum field 
operator, 

G.(>,r',«0 =E MgMWj%™ . (7 . 9) 



nlm 



ui 



Here = u> (2n + l + 3/2) — fx a describe the energy levels of the \l/ a atomic hyperfine states. 
We note that Eq. (J7.9)) does not explicitly involve the threshold energy u a or the chemical 
potential fi a , because these have been included in the effective detuning u> as defined in Eq. 
(17. 7|) . Hopefully, there will be no confusion with the label of the Fermi Matsubara frequency 
uj m = iiT(2m + 1) and the azimuthal quantum number m. 
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We substitute the diagonal Green's functions Gu in Eq. ()3.4|) and G a in Eq. (17.91) into 
Eq. (|7.8jh and carry out the Matsubara u; m -frequency summation in the usual way|46j]. After 
doing the analytic continuation to real frequencies, we obtain (at T = and uj > 0) 

I F {u) = 27r4£(2* + 1) E E \K + i,N l+j MO^i + E%- uj), (7.10) 

/ j=0 n=0 

where Ef is the energy eigenvalue given by the self-consistent solutions of the BdG equations. 

The tunneling current from molecules Im(u) = (-Tw^)) can be calculated in the same 
way. Within first order perturbation in Hj 1 , we find 

Im(u) =2t 2 M Im J drdr'(j) m (r)(t) m (r')Tl B (r,r',-uj), (7.11) 

where 

r/3 



U B (v,r',iiy n ) = - f dre^^T^^T)^,^^^')}) 

J 

Y J G 2 2(r,r',iuj m + iu n )G a (r',r,iuJ m ). (7.12) 



1 



Carrying out the a; m -frequency summation in Eq. (|7.12|) and the analytic continuation, we 
obtain 

N t N, 

I M {u) = 2ntl; ^(2/ + 1) E E K\ 2 Q(&W& + Ef t -oo), (7.13) 

l j=0 n=0 

where the matrix element is given by 



SS = E ^4+^,^^ / r 2 rfr< i (r)0 M (rX, i (r). (7.14) 



n'=0 

The analogue of these results were first worked out by Torma and co-workers [29 



54. 



55 



m, 



both for a uniform gas and a trapped gas within the LDA. Refs. [3, 0] only include the 
molecular tunneling current originating from the dissociation of excited molecules, which 
only becomes important at finite temperatures close to T c . In the case of T = 0, which we 
are considering, all the molecules are Bose-condensed. In this case, the contribution in Eq 
(17.13)) associated with the dissociation of the Bose-condensed molecules gives the dominant 
contribution to the molecular current. 

In order to illustrate the physics of the preceding expressions for Iy(u) and Im(uj), it is 
useful to consider a uniform superfluid Fermi gas at T = 0. In this case, it is convenient to 
evaluate If(u) and Im(uj) in momentum space. In a uniform gas, the two-particle Green's 

29 



function II in Eq. (|7.6|) only depends on the relative coordinate as II(r — r', —oj), so that 
Eq. ()7.6|) can be written as 

I f (oj) = 2t 2 F - Im Gn(P,iu m + iu n )G a (p,iu m ) , (7.15) 

where Gn is given in Eq. (|2.31|) . and G~ 1 (p, iuj m ) = ioj m — £ p is the single-particle Green's 
function of a free uniform Fermi gas of atoms in the hyperfme state 'a'. After doing the uj m - 
and p-summations, one finds 

I F {u>) = nt F p(Q = + - A)Q(l^^- + „). (7.16) 

K 2 OJ J U! z K 2 OJ J 

Here, p(Q) oc yT^ is normal state density-of-states given below Eq. (|6.ip . Equation ()7.16|) 
clearly has a peak at the energy gap given by oj = |A| as long as fi > (BCS region). On 
the other hand, when [i < (BEC region), the threshold energy is given by the last factor 
in Eq. (TOT)!) . 



u th = \n\ + V/i 2 + A 2 , (7.17) 

rather than = \A\. In this BEC limit, the expression in Eq. (I7.17J) approaches the 
binding energy of a molecular boson as uj^ — » 2\fi\ ~ which is twice as large as 

the threshold energy (or energy gap) of the single-particle Fermi excitations (in a uniform 
gas). As one might expect, the threshold energy of the rf-induced tunneling current If(uj) 
continuously changes from the single-particle excitation gap |A| in the BCS regime to the 
threshold of the two-particle continuum in the BEC limit. In a trap, the tunneling current 



I F (u) in LDA 



29( is given simply by the spatial integration over Eq. ()7.16|) . where A(r) 



and fi(r) = // — \vmJ^r 1 now depend on the position r. The sharp peak at the excitation 
threshold then becomes broadened, as first discussed in Ref. 2^|. 

Figure Uni shows our calculated results for the rf-induced current I F {oS) in the BCS-BEC 
crossover at T = in a trapped gas. Since we take into account the discrete energy levels 
of the harmonic trap potential (see Fig. [TUJ) . the spectrum shows rapid oscillations. In 
panel (a) descri bing the BCS region, one sees that the lowest tunneling current frequency 



is at u ~ 0.08e F 



6l|. This corresponds to the single-particle excitation gap E g discussed in 



Section VI. A broad peak is also evident, centered at uj ~ 0.3e F - This peak energy is seen 
to decrease as one enters the crossover region [panel (b) of Fig. EI]. Since the profile of the 
composite order parameter shrinks and the chemical potential \x decreases in this regime 
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(see Figs. EJandEj), the width of the combined potential well [A(r) + V^ ap (r), see also the 
inset in Fig. ITlT b)] increases. In this case, a large number of low-energy excitations appear, 
which are localized at the minimum of the potential well. These low energy excitations lead 
to the large value of If{w) at low u, as shown in Fig. Ilfif b). In contrast, once we enter the 
BEC regime, the single-particle excitation gap E g quickly becomes large (see Fig. ITTj) . This 
shows up in the rf-tunneling currents in panels (c) and (d), where the low-energy Andreev 
surface excitations no longer determine the pair threshold in If(u), even though these states 
still exist. In Fig. fTBT d). we see the threshold energy of the Fermi spectrum becomes quite 
large, at u>th — %£f, a result of being in the BEC region. This value is approximately twice 
as large as E g (see Fig. ITTj) for this value of u, and is consistent with the discussion given 
above for the case of a uniform superffuid. 

For comparison, Fig. IT7I shows the rf-tunneling spectroscopy calculated using the LDA, 
which is given by [see Eq. ()7.16jl ] 

W = < [dr P (n = i^AW! + Mr) )^ e(w - Aw)e(i^l)! + Mr)) . 

(7.18) 

Here /j(r) = // — |ma;Qr 2 + ynp(r), where ynp(r) is the Hartree term. The values of /j, A(r), 
and np(r) which are used in Eq. (|7.18j) are obtained by solving the BdG coupled equations 
in a self-consistent manner, as discussed in earlier sections. 

Comparing Fig. El with Fig. we find that when these self-consistent solutions of the 
BdG equations are used, the LDA gives a good overall approximation in the whole crossover 
regime for the rf-tunneling spectroscopy. However, it does not give the true energy gap 
(E g ~ 0.08£f *C £p) in the BCS regime, in contrast to the microscopic results shown in Figs. 
ITfiT a) and (b). Since this energy gap E g originates from the Andreev bound states, it is no 
surprise that an LDA calculation does not reproduce it. 

Although the peak energy in the tunneling current is very small in Fig. flHT b). this does 
not mean that the "average" magnitude of the composite order parameter is small. When 
we extract the contribution coming from the central region of the trap (0 < r < r c , where 
r c Rf)i If(w) has a peak at a high energy, reflecting the large magnitude of A(r ~ 0). 
This is shown in Fig. ITTT b) for the case r c = 0.3-Rf (line 1), where Ip DA (u) has a peak 
around u ~ ep- In this restricted spatial region (0 < r < 0.3-Rp); A(r) is of the order of the 
Fermi energy, as shown in the inset of Fig. ITTT b). As one increases r c , the contribution of 
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the low-energy excitations localized around the surface region of the cloud begins to "hide" 
this high energy peak. When we take r c = 0.6i?F in Fig- fT7T b) (the line 2), the peak in 
the spectrum is still dominated by the low energy excitations. It would be very useful in 
future experiments if one could measure the rf-tunneling current from the central region of 
the trap. In principle, such selective measurements could give detailed information about 
the spatial dependence of A(r). 

When the profile of the composite order parameter spreads out up to Rp (namely, when 
A(r) is large even close to the trap edge), one expects that the effect of high energy excita- 
tions with an energy gap close to A(r = 0) will dominate over low-energy excitations in the 
rf-spectrum Ip(u)). To confirm this expectation, we show a calculation in Fig. ^] using an 
ad-hoc model for A(r) shown in the inset. One finds that, as expected, the rf-tunneling spec- 
trum for this broad order parameter has a high-energy peak at u ~ A(0) ~ e-p. This kind 
of slowly decreasing A(r) might be obtained when the effective repulsive molecule-molecule 
is strong, which our calculations have ignored. We also note that our explicit calculations 
are for a narrow Feshbach resonance, in which case molecules are dominant in the crossover 
region. [Our self-consistent expression for A(r) is shown in the inset of Fig. fHT b).] In a 
broad Feshbach resonance, where Cooper-pairs are dominant in the crossover regime and the 
size of these Cooper-pairs is still fairly large, such a broad profile of A(r) may be possible. 
In the recent calculations of the rf-tunneling by Kinnunen et al.j^, a high-energy peak was 
obtained in the low temperature limit. The difference between the result obtained in Ref. 
I29I ] and ours (Figs. ITH1 and fTTj) thus seems largely due to the difference in the spatial profile 
of the composite order parameter A(r) used in the two calculations. 

As shown in Fig. and the discussion above, the rf-tunneling current If(cj) is very 
sensitive to the detailed spatial structure of A(r). This is simply because the factor oj~ 2 
in Ip{uS) [see Eq. ()7.16jl ] tends to emphasize the role of the low energy excitations of the 
superfluid gas. To extract information about the high-energy region and the magnitude of 
A(0) in the center of the trap, it is useful to consider the function If{oj) = u) 2 Iy{oj). In a 
uniform gas, when we neglect the energy dependence of the normal state density-of-states 
p(uj) for simplicity, we find that Jf(cg>) oc 0(u; — |A|), so that one can directly determine the 
magnitude of A from the energy at which the spectrum shows a sudden jump. Using the 
LDA, this discontinuity at to = | A| is broadened in a trap due to the inhomogeneity of A(r). 
However, one can still expect that If{<^) would start to decrease from an energy of the order 
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of the maximum gap, since most atoms are at the center of the trap. This behavior is clearly 
shown in Fig. O where If(u) is seen to be suppressed for to < A(r = 0). At the center 
of the trap, A(0) is of the order 1.5 ~ 2e F (see Fig. [T3J). We conclude that a plot of If(oj) 
can be used to estimate the magnitude of the composite order parameter at the center of 
the trap, even if there is a low-energy peak in If(oj). 

Figure l2*Ul shows the spectrum at T = of the molecular current Im(&) given by Eq. 
([7. 13)1 . This Bose rf-tunneling spectrum is seen to be very similar to If(&), as defined above 
(see Fig. [T9*]) . In particular, we see that the frequency where the rf-spectrum starts to be 
suppressed corresponds quite closely to the maximum value of the order parameter A(r = 0) 
at the center of the trap (as shown by the arrows in Fig. !2Uj) . In a uniform gas, one can 
calculate Im(^) explicitly, 

I m (oj) = id 2 M tf M p(n = \^-^- + A*)e(o; - A)0(i^^ + A (7.19) 

^ L OJ ' X Z OJ ' 

Thus in a uniform gas, we find that Im(oj) oc If(u), since the factor oj~ 2 in Eq. ()7.16|) is not 
present in Im(ijj). Thus, measuring the molecular dissociation current Im{w) would appear to 
give a more direct way of probing the spectrum of the high energy excitations, and hence the 



magnitude of the composite order parameter A(r = 0). In recent discussions 56(, a general 



expression for the molecular tunneling current was discussed, but it was not evaluated in 



Ref. 



This neglect was justified by the assumption that in the BCS region of interest, 



the n Umber of .olec.es was 8ma ... Fo r a broad Feshbacb _e considered in Re, Q, 
this may be correct but further studies are needed. 

Since we can eliminate the effect of low-energy excitations on the rf-spectroscopy and 
extract the information about A(0) by considering If{uj) and Im{w), it is interesting to see 
where this "hidden peak" coming from this high- energy contribution is in the spectrum of 
If(oj) in cases where it may be masked by the low energy spectral weight (as in Fig. [Ttjj) . 
For this purpose, we can simply model the If(u) spectrum using the Lorentzian form 

^ - + (7 - 20) 

The parameters C, T, and oj p can then be determined from the best fit to the calculated 
If{w) around the high-energy region u ~ A(0). As an example, such a fit is shown by the 
dashed line in Fig. The corresponding rf-tunneling spectrum is shown by the 

dashed line in Fig. [TBT b). The energy of the "hidden" broad peak in If(oj) is given by the 
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value of Up found by this procedure. The results for uj p using this procedure are plotted in 
Fig. |^ We find reasonable agreement with the experimental data for 6 Lij2?|, especially 
for the case U = 0.35ep- Our peak energy at (kpa s ) _1 = occurs at u ~ 0.3ep, which is 
in agreement with the results of the recent theoretical analysis using the LDA by Torma 



and co-workers 



291 ]. who found a broad peak at u ~ 0.3e-p near the unitarity limit. We also 
291 ] were for a broad resonance. It would be very useful 



recall that the calculation in Ref. 



to have the calculations in Ref. 29^ extended to cover the whole crossover region (results 
were only reported for v = 0.5ep). 

In the recent rf-spectroscopy data at finite temperatures of Grimm and co- workers j^. 
one finds a strong narrow peak at zero detuning as well as a broad peak at positive detuning. 
As the temperature decreases, the spectral weight shifts to the broad peak. It is argued in 
Ref. j^] that the peak at zero detuning is due to unpaired or free Fermi atoms at the edge 
of the trap, which have no energy gap. The fact that the peak is narrow is further argued to 
be evidence that these states come from the region of low density, consistent with negligible 
mean-field broadening. The pioneering theoretical work of Torma and co-workers based on 
the LDA [3] leads to the same interpretation. However, as discussed at the end of Sec. V, 
the LDA at finite temperatures incorrectly predicts a region at the edge of the trap where 
the order parameter A(r) has vanished even though the density of atoms n(r) is still finite 
(see Fig. 1 of Ref. In fact, the BdG equations show that the entire trapped gas is in 

a superfluid state below T c , with A(r) being finite everywhere where n(r) is finite. 

As we discussed in Section VI, very low energy states (with a finite but small excitation 
gap Eg) arise which are localized in the low density tail of the superfluid gas, the analogue of 



Andreev states 
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571 ] . These in-gap surface states |55j, |58| are a true signature of the Fermi 



superfluid phase, since they involve a coherent mixture of particle and hole components |5S 
In particular, these excitations exist even at T = 0. At finite temperatures below T c , 
the rf-tunneling current If(u) involves two contributions, coming from thermally excited 
Bogoliubov quasi-particles and from the Cooper-pair condensate. The former contribution 
disappears at T = 0, because all the atoms are paired and Bose-condensed, and there are 
no excitations. However, the latter contribution exists even at T = 0. 

The low-energy rf-current spectrum Ip(uj) is in fact dominated by excitations from the 
condensate associated with the low-energy Andreev or in-gap states. Indeed, we can clearly 
see the true very small excitation gap E g <C ep from the lowest energy peaks in panels (a) 
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and (b) of Fig. Torma and co-workers [29| have argued that the normal phase atoms 
at the trap edge in their LDA calculation can be viewed as a crude approximation to these 
in-gap low energy states, which arise in a more accurate theory based on the BdG equations. 
However, in an LDA analysis, the low frequency (free atom) peak is predicted to disappear 
at very low temperatures when the entire trapped gas becomes superfluid. This picture 
is different from what our microscopic calculations give, as discussed above. For the same 



reason, we feel that the data in Re 
pseudogap phase discussed in Refs. 



does not give any convincing evidence for the 

45|. 

A proper discussion of the low energy states in the edge of a trapped superfluid gas re- 
quires the kind of microscopic calculations presented in this paper. One is losing a huge 
amount of physics by thinking of the central peak (at zero detuning) simply as the contri- 
bution from unpaired atoms (or atoms feeling the effect of a "pseudogap"). This low energy 
peak at uj <C e f is due to the characteristic low energy states (~ E g <C e F ) of a trapped 
superfluid gas in the BCS region (see Fig. 110)1 . The BEC region is quite different. We recall 
from Fig. ^2 that E g rapidly increases as we pass from the BCS to the BEC regions. In the 
BEC region, the large energy gap E g is determined by the magnitude of the chemical poten- 
tial, rather than the low energy Andreev in-gap states. It would be very useful to attempt 
a higher resolution study of the "unshifted peak" in the rf-induced current measurements 



as presented in Ref. 



This could give more detailed information about the low-energy 



in-gap excitations of a trapped Fermi superfluid in the BCS region. 



VIII. CONCLUDING REMARKS 



In this paper, we have presented a detailed study of the equilibrium properties and single- 
particle excitations in the BCS-BEC crossover regime of a trapped Fermi gas with a Feshbach 
resonance, at T = 0. We extended the crossover theory developed by Leggettj^ to include 
the effect of a Feshbach resonance in a trapped Fermi gas. In our work, the composite 
order parameter A(r), the atomic density profile np(r), and the chemical potential [i are all 
calculated self-consistently. Our theory does not use the LDA, but works with the correct 
eigenstates of the harmonic trapping potential given by the Bogoliubov-de Gennes (BdG) 
coupled equations. In a uniform BCS Fermi superfluid, the single-particle excitations have 
an energy gap which is equal to the pair potential A. In contrast, in the BCS-BEC crossover 
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region, this single-particle energy gap is not directly proportional to the pair potential A. 
In a trapped Fermi superfluid, there is never a simple relation between the energy gap E g 
and the underlying spatially varying order parameter A(r). One of the themes of our paper 
is how to extract information about both these quantities. 

We showed that spatially-dependent local density and the order parameter become more 
localized at the center of the trap as one decreases the threshold energy (2v) of the Feshbach 
resonance. This reflects the fact that the character of the particles is continuously chang- 
ing, from unpaired Fermi atoms to bound states (molecules) associated with the Feshbach 
resonance. The threshold energy E g of single-particle excitations was shown to be much 
smaller than the magnitude of composite order parameter at the center of the trap in the 
crossover region, where A(r = 0) ~ Ep. This is because E g in this region is determined by 
the lowest Andreev (or in-gap) bound states j^] near the bottom of the combined potential 
well composed of the off-diagonal pair potential and the trap potential. We have emphasized 
that these states are a characteristic signature of a trapped superfluid Fermi gas and hence 
are of special interest. 

We also used our results for the single-particle excitation spectrum to discuss recent rf- 
tunneling experiments. As discussed recently 0, 13 the data at finite temperatures in the 
crossover region can be usefully described in terms of a narrow unshifted peak (u ~ 0) and 
a broad peak at a detuning frequency u comparable to the expected pair potential A(r = 0) 
at the center of the trap. While LDA calculations for the case v = O.5£f0] appear to 
confirm this kind of rf-spectrum, our present calculations based on the explicit solutions 
of the BdG equations at T = lead to somewhat different predictions about the low uj 
region. We have verified that our results are essentially reproduced by an LDA calculation 
if we base it our self-consistent values of rip(r) and A(r) given by the BdG equations. The 
major difference between our results and those in Refs. 0, 13 is that we find a strong 
low frequency contribution to the fermionic tunneling current in the BCS-BEC crossover 
region. However, we find that the rf-tunneling current If(uj) is very dependent on the 
precise spatial dependence of A(r). This dependence is good news since it means that fits 
to the rf-tunneling data may be used in the future to obtain detailed information about the 
spatial dependence of the composite order parameter A(r) in the crossover region. 

Further work is needed to clarify the role of the Andreev bound states in determining the 
low frequency peak in the rf-spectrum data. These single-particle states continue to exist 
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even at T = 0, and are clearly left out of an LDA type calculation, which is based on the 
results for a uniform Fermi superffuid. 

In Sec. VII, we showed that the low frequency peak in the rf-spectrum was due to 
contribution from the edge of the trap, where A(r) is very small [see Fig. IT7T b)]. Subtracting 
out this "surface" contribution, the remaining rf-spectrum was peaked at high energies 
comparable to the value of A(r = 0) at the center of the trap. These peak energies were 
in reasonable agreement with the observed position of the broad peak in the tunneling data 
reported in Ref. 27[, as shown in Fig. 

We have also evaluated the rf-spectrum due to the current produced by dissociation of 
molecules (see Fig. I2UJ) and found that it did not have an unshifted component at u ~ 0. 
Thus it gives a more direct measurement of A(r = 0) at the center of the trap. This 
contribution, which was discussed but not explicitly evaluated in Ref. deserves further 
study, both experimentally and theoretically. In future work, we hope to extend our present 
calculations to finite temperatures. 

Besides the sing le-particle properties, collective excitations in trapped Fermi 
superfluidsj^, 0, Q] are also of great interest in the crossover region. The single-particle 
Green's function evaluated in the present paper form the basis for the calculation of the 
collisionless collective modes in a trapped Fermi gas using linear response theory. In a 
future paperQ, we wil. exteud the approach giveu iu Ref. 3 for a uuiform gas aud dis- 
cuss the quadrupole and monopole modes in the BCS-BEC crossover region [3 lj. A detailed 
discussion of the Kohn mode has been recently given in Ref. j^. 

As discussed in Sec. VII, rf-tunneling spectroscopy experiments can give information 
about the true quasi-particle excitation gap E g . However, very high resolution would be 
necessary in the BCS and crossover regime because E g ~ u <C £f (see Figs. ITTIandlTBj). As a 
result, an interesting problem remains as to how to measure E g in this region. An alternative 
method might be through the study of the collective mode frequencies, which typically have 
a very low energy comparable to the trap frequency ujq and are, in fact, bounded by the 
two-particle continuum at 2E g . We show a plot of the calculated frequency of the monopole 
mode in Fig. E21(the details are discussed in Refs. j^jH^I). The monopole mode frequency 
is seen to be suppressed, so that it always lies below the two-particle continuum at 2E g . In 
particular, the suppression is quite striking around the region where E g shows a minimum as 
a function of v. This effect of the two-particle continuum on the monopole mode frequency 
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appears to be an attractive way of obtaining information about the single-particle excitation 
gap Eg in a trapped superffuid Fermi gas in the crossover region. 

Acknowledgments 

Y. O. was financially supported by a Grant-in- Aid for Scientific research from the Ministry 
of Education, Culture, Sports, Science and Technology of Japan, as well as by a University 
of Tsukuba Research Project. A. G. acknowledges a research grant from NSERC of Canada. 



38 



[1] Y. Ohashi and A. Griffin, Phys. Rev. Lett. 89, 130402 (2002). 
[2] Y. Ohashi and A. Griffin, Phys. Rev. A 67, 033603 (2003). 
[3] Y. Ohashi and A. Griffin, Phys. Rev. A 67, 063612 (2003). 

[4] J. N. Milstein, S. J. J. M. F. Kokkelmans and M. Holland, Phys. Rev. A 66, 043604 (2002). 
[5] A. Perali, P. Pieri, and G. C. Striati, Phys. Rev. A 68, 031601 (2003). 

[6] E. Timmermans, K. Furuya, P. W. Milonni and A. K. Kerman, Phys. Lett. A 285, 228 (2001). 
[7] M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo and R. Walser, Phys. Rev. Lett. 87, 

120406 (2001). 
[8] D. M. Eagles, Phys. Rev. 186, 456 (1969). 

[9] A. J. Leggett, in Modern Trend in the Theory of Condensed Matter edited by A. Pekalski and 

J. Przystawa (Springer Verlag, Berlin, 1980), p. 14. 
[10] P. Nozieres and S. Schmitt-Rink, J. Low. Temp. Phys. 59, 195 (1985). 
[11] A. Tokumitsu, K. Miyake and K. Yamada, Phys. Rev. B 47, 11988 (1993). 
[12] M. Randeria, in Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke and S. 

Stringari (Cambridge University Press, N.Y., 1995), p. 355. 
[13] C. A. R. Sa de Melo, M. Randeria and J. R. Engelbrecht, Phys. Rev. Lett. 71, 3202 (1993). 
[14] J. R. Engelbrecht, M. Randeria, and C. A. R. Sa de Melo, Phys. Rev. B 55, 15153 (1997). 
[15] R. Haussmann, Self- consistent Quantum-field Theory and Bosonization for Strongly Correlated 

Electron Systems (Springer- Verlag, Berlin, 1999) Chap. 3. 
[16] T. Loftus, C. A. Regal, C. Tichnor, J. L. Bohn, and D. S. Jin, Phys. Rev. Lett. 88, 173201 

(2002). 

[17] K. Dieckmann, C. A. Stan, S. Gupta, Z. Hadzibabic, C. H. Schunck, and W. Ketterle, Phys. 

Rev. Lett. 89, 203201 (2002). 
[18] C. A. Regal, C. Tichnor, J. L. Bohn, and D. S. Jin, Nature 424 (2003) 47. 
[19] K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. Rev. Lett. 91 (2003) 080406. 
[20] M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 (2003). 

[21] S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, S. Riedl, C. Chin, J. Hecker Denschlag, 

and R. Grimm, Science, 302, 2101 (2003). 
[22] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004). 



39 



[23] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J. Denschlag, and R. Grimm, 

Phys. Rev. Lett. 92, 120401 (2004). 
[24] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Raupach, A. J. Kerman, and W. 

Ketterle, Phys. Rev. Lett. 92, 120403 (2004). 
[25] J. Kinast, S. Hemmer, M. Gehm, A. Turlapov, and J. Thomas, Phys. Rev. Lett. 92, 150402 

(2004). 

[26] T. Bourdel, L. Khaykovich, K. Cubizolles, J. Zhang, F. Chevy, M. Teichmann, L. Tarruell, S. 

J. J. M. F. Kokkelmans, and C. Salomon, Phys. Rev. Lett. 91, 020402 (2003). 
[27] C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, J. Denschlag, and R. Grimm, 

Science 305, 1128 (2004). 
[28] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. Chin, J. Denschlag, and R. Grimm, 

Phys. Rev. Lett. 92, 203201 (2004). 
[29] J. Kinnunen, M. Rodriguez, and P. Torma, Science 305, 1131 (2004). 



[30] Y. Ohashi and A. Griffin, IcomTmat /0402031 



[31] Y. Ohashi, in Workshop on Ultracold Fermi Gases, Levico (Trento), Italy, March 4-6, 2004. 

This paper is available as a pdf-file at: |http://bec.science.unitn.it /fermi04 
[32] T. Kostyrko and J. Ranninger, Phys. Rev. B 54, 13105 (1996), and references therein. 
[33] P. Pieri, and G. Strinati, Phys. Rev. B 61, 15370 (2000). 

[34] D. Petrov, C. Salomon, and G. Shlyapnikov, Phys. Rev. Lett. 93, 090404 (2004). 

[35] A. Griffin, Excitations in a Bose-Condensed Liquid (Cambridge University Press, Cambridge, 

1993), Chapters 3 and 5. 
[36] Y. Ohashi and A. Griffin, in preparation; see also Ref. [^|. 
[37] G. Bruun, Y. Castin, R. Dum, and K. Burnett, Eur. Phys. J. D 7, 433 (1999) 
[38] G. Bruun, and B. R. Mottelson, Phys. Rev. Lett. 87, 270403 (2001). 
[39] G. Bruun, Phys. Rev. Lett. 89, 263002 (2002). 
[40] G. Bruun, and H. Heiselberg, Phys. Rev. A 65, 053407 (2002). 

[41] P. G. de Gennes, Superconductivity of Metals and Alloys (Addison- Wesley, NY., 1966), Chap- 
ter 5; T. Tsuneto, Superconductivity and Superfluidity (Cambridge University Press, 1998), 
Chapter 3. 

[42] Since H^p^ is independent of the quantum number m, we need only consider H$ FB . 

[43] J. Stajic, J. Milstein, Q. Chen, M. Chiofalo, M. Holland, and K. Levin, Phys. Rev. A 69, 



40 



063610 (2004). 

[44] A. Perali, P. Pieri, L. Pisani, and G. Strinati, Phys. Rev. Lett. 92, 220404 (2004). 
[45] For a review, see Q. Chen, J. Stajic, S. Tan, and K. Levin, cond- mat/0404274 



[46] G. D. Mahan, Many Particle Physics, Second Edition (Plenum Press, N.Y., 1993), Chapter 3. 
[47] J. R. Schrieffer, Theory of Superconductivity (Addison- Wesley, Redwood, California, 1964), 
Chapter 7. 

[48] See, for example, A. Fetter and J. Walecka, Quantum Theory of Many-Particle Systems 

(McGraw-Hill, N.Y., 1971), p.l28ff. 
[49] We note that a« is different from a many-body scattering length related to a many-body t- 

matrix T(q, u;)[j| in the low-energy limit, defined by Airaf^ /m = T(q = 0, lo = 0). This always 

diverges below T c due to the presence of the Goldstone mode at zero frequency 3j. 
[50] The slight difference between the two is due to the small difference of the chemical potential 

determined by the equation for the total number of atoms, as shown in Fig. Ufa). 
[51] M. L. Chiofalo, S. J. J. M. F. Kokkelmans, J. N. Milstein, and M. J. Holland, Phys. Rev. Lett. 

88, 090402 (2002). 
[52] P. G. de Gennes and D. Saint James, Phys. Lett. 4, 151 (1963). 

[53] This "excitation" energ y g ap is measured relative to the value of the chemical potential /i. In 

contrast, Fig. 3 of Ref. |30j gives a plot of {i + E g . 
[54] P. Torma, and P. Zoller, Phys. Rev. Lett. 85, 487 (2000). 

[55] G. Bruun, P. Torma, M. Rodriguez, and P. Zoller, Phys. Rev. A 64, 033609 (2001). 
[56] J. Kinnunen, M. Rodriguez, and P. Torma, Phys. Rev. Lett. 92, 230403 (2004). 
[57] M. Baranov, JETP Lett. 72, 385 (2000). 
[58] H. Heiselberg, |cond-mat/0406714[ 

[59] D. Luxat, and A. Griffin, Phys. Rev. A 65, 043618 (2002). 

[60] C. Choi, Y. Japha, and K. Burnett, Phys. Rev. A 61, 063606 (2000). 

[61] This value of the excitation threshold is rather larger than E g ~ 0.04ep given in Fig. 111! 
The difference is simply due to the fact that the threshold energy (= E^) in the tunneling 
current If{u) is determined by E t ^ = Mm[£^ + E^] = Min[^] + E g: under the condition 
that energy of the ^ a -state [^ = h~u!o(2n + I + 3/2) — n] is positive. In panel (a) of Fig. 
1161 the chemical potential is /i = 27.56^0- Choosing the state with quantum numbers (n, I) 
such that 2n + I = 27, Mm[£^] ~ 0.94cj ^ 0.03e F . This gives a tunneling threshold energy 

41 



E th = E g + Min[^\ ~ 0.04e F + 0.03e F = 0.07e F ^ 0.08e F in agreement with Fig.EEfa). 



[62] H. Heiselberg, cond-mat/0409077 ; references to the literature are given here. 



[63] Y. Ohashi, cond-mat/0407379 ; Phys. Rev. A, in press. 



42 



FIG. 1: (a) Superfluid transition temperature T c as a function of the scattering length (/c^a s ) _1 
for a uniform gas. The solid line shows the case of a narrow Feshbach resonance (g r \/n ^ £f)j 
while the dashed line is for the case of a broad Feshbach resonance (g T y/n S> £p), where n is the 
he peak structure around (/cfOs) -1 ~ 0.25 is an artifact of the approximation 



density of atoms. 

n 

and is not intrinsic [15J]. The solid circle shows T c as a function of the two-particle scattering length 
a^ b in the case of g v yjn = 20£f (broad Feshbach resonance). We take the cutoff oj c = 2ep- (b) 
The number of Fermi atoms N-p and the number of Bose molecules Am at T c in the BCS-BEC 
crossover. (Nm includes both Cooper-pairs and real two-body molecules.) The total number of 
particles is given by N = Np + 2 Am- Irrespective of whether one is dealing with a broad or narrow 
Feshbach resonance, one notes that bound states (Nm) become dominant when (fepOs) 1 <; 0. 

7 IG. 2: (a) Chemical potential in the BCS-BEC crossover for a uniform gas at T c , based on Ref. 
jj. The parameters are the same as in Fig. ^ Panels (b) and (c) show detailed character of 
particles at T c for a uniform gas in the cases of a narrow Feshbach resonance and broad Feshbach 
resonance, respectively. N^ >0 is the number of Feshbach molecules with a finite lifetime, and N^~ 
is the number of stable Feshbach molecules. The number of stable Cooper-pairs is N%- ; JVg= is the 
contribution from the scattering states or Cooper-pairs with a finite lifetime. For precise definitions 
of iV B and N c , see Ref. Q. 



FIG. 3: (a) Superfluid transition temperature T c as a function of the scattering length (kpas) -1 
for a trapped gas. The cases of a narrow Feshbach resonance and broad Feshbach resonance are 
respectively shown by the solid line and dashed line, respectively. In this figure, n = N/R F , and we 
take the cutoff lo c = 2ef- (b) The number of Fermi atoms A^f and the number of Bose molecules 
A^m at T c in a harmonic trap. 
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FIG. 4: (a) Comparison between the bare attractive interaction — U e s between Fermi atoms in the 
gap equation [as defined in Eq. (|2.1U[) ] and the interaction Ana s /m defined in Eq. (|4.7|) . We take 
the cutoff frequency lo c = 161.5a;o and the Fermi energy ep = 31.5u;o i n this and later figures, (b) 
Relation between the scattering length a s and the bare threshold energy 2v, for the parameters we 
use in this paper. The Feshbach coupling strength is taken as g T = 0.2ep- 

FIG. 5: Density profile in the BCS-BEC crossover. np(r) and riM(r) represent the Fermi atom 
density and Bose molecule density, respectively. n(r) = np(r) + 2nyi(r) is the total density profile, 
r is normalized by Thomas-Fermi radius Ftp = J 2ep / moj\ for a free Fermi gas in a trap. We take 
g T = 0.2ep, U = 0.35£f and ojqb = <^o- These values are are also used in Figs. Eld and|Hl 

FIG. 6: Fermi chemical potential [i in the BCS-BEC crossover at T = 0. The inset shows fi as a 
function of the threshold energy 2v. The dashed line in the inset is the lowest energy of molecular 
excitation spectrum in a trap. 

FIG. 7: Spatially-dependent composite order parameter A(r) in a trap in the BCS-BEC crossover. 
The Cooper -pair order parameter A(r) and molecular condensate (fi m (r) are also shown. 

FIG. 8: Expansion coefficient a n of the molecular condensate wavefunction 4>m ( r ) = 
-7= J2 n a nRno{ r ) given by Eqs. (|2.21|) and (|2.22|) . BEC in a non-interacting Bose gas corresponds 
to all a n = for 

FIG. 9: Comparison of the calculated atomic density profile n(r) and composite order parameter 
A(r). The results are normalized by the values at the center of the trap, (a) BCS regime, and (b) 
BEC regime. The inset compares n(r) with A(r) 2 in the BEC regime. 

FIG. 10: Single-particle density of states (DOS) at T = 0. We take g T = 0.2e F and U = 0.52e F . In 
calculating the DOS, we have introduced a small imaginary part (r = 0.005u;o) in the eigenenergies. 
The fine structure in DOS is due to the discrete levels of harmonic trapping potential. The peaks 
are the single-particle Bogoliubov quasi-particle energies in a trapped gas. Note that fi < in 
panel (c). 
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FIG. 11: Single-particle excitation gap E g appearing in the density of states N(u>) in the BCS-BEC 
crossover in a trap. The dashed line shows in the negative fj, region [(/jf^s) 1 > 0.65]. The 
region around the minimum of E g shown in the inset is discussed in the text. We take g r = 0.2ep 
and U = 0.52e F . 

FIG. 12: A plot of the calculated number of atoms (Np) and molecules (Nm) at T = 0. The total 
number of atoms is given by N = Np + 2JV M . We find iV M ~ N/2 = 10912/2 in the negative (i 
region (see also Fig. |f)J). We take g T = 0.2ep and U = 0.52sf- 

FIG. 13: Change in the maximum value of composite order parameter A(0) at the center of the 
trap in the crossover region. 

FIG. 14: Localized density of states N(u, r) given by Eq. (|3,9|) in the crossover regime. The peaks 
correspond to the Andreev or in-gap states. The inset in panel (b) shows the combined potential 
well consisting of the composite "off-diagonal" pair potential A(r) and "diagonal" trap potential 
measured from the chemical potential, given by V^ ap (r) = (V t f &p (r) — fJ>)@(V t f ap (r) — fx). Actually, 
atoms also feel the Hartree potential —^np(r), but this is not plotted in the inset. 

FIG. 15: (a) Combined potential well formed by the off-diagonal pair-potential A(r) and the 



diagonal trap potential V t f ap (r). The combined potential is given by the sum of A(r) and V e g 



trap 



r = 



(Vtrap( r ) — fi)@(V t f ap — (j,). The low-energy states appear around the bottom of this effective 



potential wel 
bound states 



52, 



Simplified model of the combined potential well which gives rise to Andreev 



FIG. 16: The rf-induced current If(w) in the BCS-BEC crossover region at T = 0. We take 
g T = 0.2ep and U = 0.52ep, and introduce a finite imaginary part T = 0.5ujq to the Bogoliubov 
eigenstate energies to smooth out the results. The small but finite intensity at u = in panels 
(a) and (b) is due to this imaginary part. The dashed line in panel (b) is obtained by fitting Eq. 
(j7.2Uf) to If(uj) = u) 2 Ip{uj) around u> ~ A(r = 0) (see Fig. IT§)) to extract the high-energy "hidden" 
peak associated with the large composite order parameter A(r = 0) at the center of the trap. The 
rapid oscillations in the tunneling spectrum, also seen in Figs. 1191 and QUI originates from discrete 
quasiparticle energy levels in a harmonic trap potential. 
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FIG. 17: The rf-induced current If(w) evaluated using the LDA. The parameters are the same 
as Fig. ED In calculating the spectrum, we have used the correct values of np(r), A(r), and fi 
obtained from the solutions of the BdG coupled equations. In panel (b), (1) shows the spectral 
contribution from the trap region < r < 0.3Rp, and (2) shows the contribution from the larger 
trap region < r < 0.6Rp. These results clearly show that tunneling from the low-energy surface 
states is the origin of the large low frequency peak. The inset in panel (b) shows the profile of the 
computed composite order parameter A(r) used for v = 0.6£f- 

FIG. 18: The solid line shows the rf-induced current If{u) based on the ad-hoc "broad" composite 
order parameter A(r) shown in the inset. This has the same maximum value A(0) at the center 
of the trap as the correct order parameter shown in the inset of Fig. I17f b). However, the width is 
broader. In this calculation, we have also used an ad-hoc broad density profile n-p(r), with a width 
of the same order as A(r). The dashed line shows the rf-spectrum plotted in Fig. I17f b). based on 
the self-consistent values of A(r) and np(r). 

FIG. 19: Spectrum of If(lo) = ll> 2 If(uj) versus the detuning frequency u, for v = 0.6£f- Parameters 
are the same as in Fig. I16f b). The arrow shows the value of A(r = 0) at the center of the trap. 
The dashed line shows a fit to Eq. (|7.20j) around the high energy region (uj ~ A(0)). 

FIG. 20: Molecular dissociation current Im(u) as a function of the effective detuning frequency 
uj. The arrow shows the value of the composite order parameter A(0) at the center of the trap. 
Parameters are the same as in Fig. 1191 

FIG. 21: Peak energy in the rf-tunneling current. In the crossover regime (0.2 < /u/ef < 1) ; the 
current is peaked at an energy evaluated using Eq. (|7.20f> . The solid line shows A(0) in the case 
of U = 0.35£f (see also Fig. I13|) . Experimental data (open circles) are taken from Figs. 1 and 2 of 

Ref. bai. 
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FIG. 22: Calculated frequency of the monopole mode (L = and n = 1) in the BCS-BEC 
crossover region at T = 0. In this figure, 2E g describes the threshold of the two-particle continuum 
spectrum. The monopole mode is suppressed below this two-particle continuum. The collective 
mode frequency is obtained from the pole of density-density correction function calculated in the 



generalized random phase approximation 
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36]. 
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Fig. 10: Ohashi and Griffin 
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Fig. 16 Otoastoi and Griffin 
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